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Abstract 

Even distribution of carriers allows to maximize 
optical gain of the Optically-Pumped Vertical Exter- 
nal Cavity Surface Emitting Laser. In this paper we 
show how to distribute the quantum wells and block- 
ing layers in order to compensate the exponential de- 
cay of the pumping beam intensity. Our model says 
whether it is possible at all (for an assumed length of 
the device) and, if it is, allows to find positions of the 
blocking layers. No iterations nor numerical calcula- 
tions more sophisticated than a standard calculator 
can do are required to use the model. 



1. Introduction 

The Optically Pumped Vertical External Cavity 
Surface Emitting Lasers (OP-VECSELs) are able 
to emit high quality beams of multi-watt pow- 
ers [T], thus combining the most important advan- 
tages of Edge-Emitting and Vertical Cavity Sur- 
face Emitting Lasers (EELs and VCSELs). More- 
over their external cavity may contain additional 
elements like non-linear crystals or semiconductor 
saturable absorber mirrors (SESAMs), which may 
be used for instance for frequency doubling and 
short pulse generation. Usually, optical gain is pro- 
vided by several quantum-well (possibly double- or 
multi-quantum-well) active regions, located at the 
successive anti-nodes of the standing wave. The cav- 
ity is formed by a Distributed Bragg Reflector (DBR) 
and an external mirror [2J. 

The pumping beam penetrates the device in the 
direction nearly perpendicular to the layers. Absorp- 
tion of the pumping beam generates carriers neces- 
sary to achieve optical gain, but also causes expo- 
nential decay of the intensity in the deeper regions. 
If we simply placed, at each anti-node, the same 
number of wells, the gains provided by the active re- 
gions would differ significantly. Since optical gain is 
a concave function of carrier concentration (roughly 
~ log(A/"/A/"o), where J\f is the carrier concentration 
and Afo is the transparency concentration), the high- 
est possible total gain (for a fixed, arbitrary number 
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of the carriers) is highest when all the wells provide 
the same material gain. If the temperatures of the 
wells can be assumed to be equal, we should try to 
have the same carrier concentration in each well. The 
temperature rise in the VECSEL can be high, most 
of the temperature drop takes place in the substrate 
and DBR as these part are much thicker than the 
active part. Generally, we want to direct the same 
number of carriers to all the wells. 

In order to do so, we have to increase the vol- 
ume from which the distant wells collect the carriers 
or increase the number of the wells in the stronger 
pumped regions. To do the first thing one can use 
so called blocking layers — thin wide-gap layers which 
block carrier diffusion (in case of AlAs blocking layers 
in GaAs, the thickness of a few nanometers is suffi- 
cient to block the carrier diffusion). They define the 
volume from which the well(s) between them collect 
the carriers. [3J. The question how to place them is 
not trivial, mainly because positions of the wells are 
restricted to the anti-nodes, which strongly restricts 
positions of the blocking layers. 

If there are two or more anti-nodes between two 
subsequent blocking layers, one has to take into ac- 
count the carrier diffusion in order to find the actual 
carrier concentration in the two (or more) active re- 
gions. This makes the analysis more complicated [3J. 
Our goal is to build an analytical model which al- 
lows to construct the desired scheme without using 
complicated calculation. 



2. The model 

As we mentioned, in order to avoid consideration 
of carrier diffusion, we restrict our interest to the 
case, where in each segment (area bounded by the 
neighbouring blocking layers) there is only one active 
region (we treat the DBR as the Oth blocking layer, 
see Fig. [TJ, and the active regions are placed at each 
anti-node. The window layer acts as the last blocking 
layer and due to the optical reasons must be located 
at an anti-node. This means that the last segment 
must be thicker than d (see Fig. [I]). Therefore, in 
order to have the same number of generated carrier 
per a QW, we have to put more wells in this segment 
and treat it in a different manner in our calculations. 
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In order to obtain a handy result we base our model 
on the following simplifying assumptions: 

1. Widths of the wells and blocking layers are negli- 
gible. More precisely, their presence (the blocking 
layers do not absorb the pump, on the other hand 
the wells have higher absorption than the adjacent 
bulk material) introduces a tolerable error. 

2. The blocking layers block totally carrier diffusion 

3. We neglect reflection of the pumping light from 
the DBR and from the blocking layers 

4. Carrier losses in the absorbing barrier are negligi- 
ble (relative to the losses in the QWs) 

In Fig. Q] a scheme of the VECSEL is presented. 
Distance d must be a multiple of a half of the emitted 
wavelength. Usually dm X/ (2n r ), where A is vacuum 
wavelength and n r is refractive index of the barrier. 
Small deviations from the exact equality come from 
the presence of the wells and blocking layers. In our 
scheme the pumping light comes from the right, so 
under our assumptions the pumping wave intensity 
is described by the following formula: 

I{z) — Iq exp(az) (1) 

where a is the absorption in the barriers, Jo is a nor- 
malisation constant, defining the pump power. Num- 
ber of the carriers generated in n-th segment is simply 

P n = I (exp(az n ) - exp(az n _i)) (2) 

We assume that in segments 1,2, ... ,N there are the 
same number of wells (in this paper — just one well) , 
and in the last segment, N + 1, we put k wells, we 
want to distribute the blocking layers such that 

Pi = Pi = ■ ■ ■ = Pn = \ p N+i (3) 

As we have k times more wells in the last segment, we 
want to have k times more carriers generated there. 
Possible values of numbers N and k will be deter- 
mined in our calculations. Position z n can be written 
as (see Fig. [1]) 

z n — nd + S n (4) 

Then we have 

P 1 =I [exp(a(d + S 1 ))-l] = 

= Iq [ exp(ad) exp(a<5i) — l] 
P 2 = I [exp(a(2d + S 2 )) - [exp (a(d + S±))] = 

= Iq [exp(2ad) exp(aS 2 ) — exp(ad) exp(a<5i)] 
P 3 = I [cxp(a(3d + 5 3 )) - [exp(a(2d + S 2 ))] = 

= Jo [exp(3ad) exp(a5 3 ) — exp(2ad) exp(a<52)] 

Let us introduce the following symbols: 

a = exp(-ad) x n — exp(a<5„) (5) 



Note that < a < 1, regardless of a and d. Using 
these symbols we can write the system of equations 
<J3]) as: 

= x 2 — 2ax\ + a 2 
= Xs — 2ax 2 + a 2 x\ 
= X4 — 2ax 3 + a 2 x 2 

(6) 

= xjy — 2axN-\ + cl 2 xn- 2 

= xn + i — (k + 1)cixn + ka 2 XN^i 

— = xn+i because 5n + i = d 
a 

As one can see, thanks to the extraordinary prop- 
erties of exponential function we got something as 
simple as a system of N + 1 linear equations with 
N + 1 unknowns. Because we do not know what are 
the values of N and k, it is convenient to consider first 
only equations concerning first N segments, i.e. those 
defined by the blocking layers which position we can 
choose. Thus we consider the system of N — 1 equa- 
tion with N unknowns: 

x 2 — 2ax\ + a 2 = 
X3 — 2ax 2 + a 2 xi = 
X4 — 2ax 3 + a 2 x 2 — 

'■ (7) 

x _/v — 2clxn-i + a 2 x.N- 2 = 

If we treat one of the unknowns as a parameter we 
can solve the system. As the parameter we choose 
x\ as the first blocking layer must be always present. 
This way we obtain: 

x 2 = a(2xi — a) 
x 3 — a 2 (3xi — 2a) 

'■ (8) 
x N = a N - 1 (Nxi - (N — l)a) 



Although the system ([7]) has always the solution in 
real numbers, we must check if the solution fulfils the 
additional conditions, i.e.: 

Kx n <l/a Vn= 1,2,..., N (9) 

The above conditions say simply that < 8 n < d. 
It assures that in each segment there is one active 
region. 
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Figure 1. Scheme of wells and blocking layers distribution in a VECSEL. Wells are blue (downwards), blocking 
layers red (upwards). Numbers 1, 2, . . . ,N + 1 denote the segments. 



Substituting © to ® we get: 



1 < Xl < - 

K a+ «) <xi< K a+ ^ 
U 2a+ £> <xi< U 2a+ ^ 



: (10) 



iV 



If the above inequalities are inconsistent, it is impos- 
sible to build a system of N active regions (of the 
same number of QWs), having equal carrier concen- 
trations. 

The above system can be significantly simplified. 
Let us denote 



L n — 



K (n - i)a+ ^) (u) 

(12) 



1 I , x 1 

Rn = -[{n- la + - 

n \ a* 



being just the left and right hand side of the n-th 
inequality. Basic calculations show that: 

L n +i - L n = * (a n+1 - (n + l)a + n) 
n(n + lja" 

(13) 

The sign of the difference is determined by polyno- 
mial l n (a) — a n+1 — (n + l)a + n. It is easy to sec 
that l n (l) = 0, l' n {a) = (n+l)(a"-l) < OVa g [0,1]. 
It means that l n {a) < Va £ [0,1] (in our case 
< a < 1), and hence 



L x < L 2 < ■ ■ ■ < L 



N 



(14) 



Thanks to that property we can replace all the in- 
equalities L n < xi, n = 1, 2, . . . , N by one: Ln < x\. 
There is no similar properties of the right hand sides. 



Now we can write the system of inequalities (fT0"| in 
a more compact way: 



1 

N 



(N - l)a + 



1 



,JV-1 



< Xl < 



mm 

n=l,...,JV 



(n - l)a + 



K{N) (15) 



Which of the right hand sides is the actual minimum 
depends on a (i.e. on the material absorption and 
distance between the active regions). 

If for a certain N Ljv > TZ(N), it is impossible 
to build N (and any greater number, since Ln is in- 
creasing and 1Z(N) is non-increasing with N) equally 
pumped segments. This way we can find the possible 
values of N. When we choose one of them, we can 
return to the full system of equations ^ . Solving it 
we can find values of all x%, X2, ■ ■ ■ , a; at, but now we 
are interested only in x\: 



1 



Xl 



(N + k)a N + 1 



N + k-1 

1 N + k 



(16) 



If, for the chosen N, we can find a natural number k 
such that xi calculated with the above formula fulfils 
condition: 

L N < xi < K(N) (17) 

we know that we can build N one-well segments 
concluded by one fc-well segment, having equally 
pumped quantum wells. If there is no such k, or 
it is too high from the practical point of view, one 
has to decrease number N and repeat the procedure. 

Finally, when we have suitable N, k and hence xi, 
we can calculate all the other X2, X3, . . . , xjv using 
and then the actual positions z n ,n — 1, 2, . . . , N + 1 
of the blocking layers using formulae (0]) and ([5]). 



3. Example 

Let us consider an important example: a VECSEL 
with GaAs barirers, emitting at 980 nm, pumped by 
a 808 nm laser. Assuming a = 1.3 • 10 4 1/cm 0], and 
n r = 3.52 we get the following parameters: 



139.2 nm 



a « 0.8345 



(18) 
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Of course in a real calculations one should use more 
precise values, because the optical properties of the 
structure are more sensitive to the distances between 
the wells. First we will check whether it is possi- 
ble to have seven one- well segments. We calculate, 
using (fTTj): 

L-i = 1.13837 72,(7) = min R n = R 3 = 1.12997 

n— 1,....7 

(19) 

We see that L 7 > 72.(7), so we cannot construct seven 
such segments. But since 

L 6 = 1.10730 < 11(6) = R 3 = 1.12997 (20) 

we can find out if we can close such 6-segment se- 
quence by the multi-well one. Using (fTB")) for N = 6 
we calculate: 

xx (k = 3) = 1.13612 > 72(6) (21) 
xi(k = 4) = 1.10595 < L 6 (22) 

Since x\ is a decreasing function of fc, there is no pos- 
sibility to find a suitable k. It means that we cannot 
build the final segment with an integer number of 
wells such that it has the same carrier concentration 
as in all the others. The non-integer values of k could 
have physical meaning as A; is actually the ratio be- 
tween number of wells in the last segment and num- 
ber of wells in the other ones. If we assumed the other 
segments to contain not one, but two wells each, we 
could consider k = 3.5. But it would mean that we 
have to put as many as seven quantum wells in the 
last segment, which is too many, because the periph- 
eral wells would be placed far from the anti-node of 
the standing wave. 

If not 6, let us try N = 5. Now we have L$ = 
1.08005, 72(5) = 1.12997. For the new N: 

Xl (k = 2) = 1.13837 > 72.(5) (23) 
L 5 < xi(fc = 3) = 1.10038 < 72.(5) (24) 
x x {k = 4) = 1.07083 <L 5 (25) 

Number x\ (k = 3) fulfils all the conditions, so we can 
build 5 one- well segments and one 3- well segment on 
the top, with all the 8 wells equally pumped. Now we 
calculate X2, ■ ■ ■ ,xn using ([S]), extract S n given by: 

S n = ^ (26) 
a 

Finally we get z\, . . . , z/v+i from formula Posi- 
tions of the wells (which are independent on X\) and 
blocking layers are presented in table [TJ One should 
remember that positions of the wells (nd) and hence 
numbers z n are only approximations. The reliable 
values are S n — positions of the blocking layers rela- 
tive to the adjacent quantum well (see Fig. [I]). 

Looking at the values in the table one can see 
that the lowest distance between a well and a block- 
ing layer is 37 nm. This is a safe distance from the 



Table 1. Positions of the wells and blocking layers in 
the 6-segment scheme, rounded to full nanometers 





[nm] 


n 




S n 


nd 




(n + l)d 


1 


1.10038 


74 


139 


213 


278 


2 


1.14012 


101 


278 


379 


418 


3 


1.13656 


98 


418 


516 


557 


4 


1.10293 


75 


557 


632 


696 


5 


1.10578 


37 


696 


733 


835 


6 


1.19837 


139 


835 


974 


974 



technological point of view, even taking into account 
non-zero thicknesses of the wells and the blocking 
layers. The total thickness of the absorbing area is 
7d, exp(— Ida) « 0.28, and hence over 70% of the 
pumping power is absorbed. As the blocking layers 
are generally located near the nodes of the stand- 
ing wave, and their thickness can be as low as a few 
nanometers, their presence modifies the optical prop- 
erties of the resonator in a very limited degree. 



4. Summary 

We have shown an efficient and simple way to de- 
sign the VECSEL structure such that the carrier con- 
centration, and hence optical gain (with the assump- 
tion that the temperature differences between the 
wells are not high), is equal in all the quantum wells. 
This configuration gives the highest modal gain for 
given number of carriers, so is highly desirable. 

We presented an example of a GaAs-based struc- 
tures with 6 active regions in subsequent anti-nodes 
of the standing wave. This is the longest possible 
design in which one use 1-well segments except the 
last one. Our calculations can be easily modified to 
describe segments with different number of wells. In 
this case longer absorbing areas can be achieved, if 
necessary. However, the longer the area is, the higher 
temperature differences appear there, which spoils 
the desired gain uniformity. 
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